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ABSTRACT 

The formation and evolution of the circumstellar disk in the collapsing molec- 
ular cloud is investigated from the prestellar stage resolving both the molecular 
cloud core and the protostar itself. In the collapsing cloud core, the first (adi- 
abatic) core appears prior to the protostar formation. Reflecting the thermo- 
dynamics of the collapsing gas, the first core is much more massive than the 
protostar. When the molecular cloud has no angular momentum, the first core 
falls onto the protostar and disappears several years after the protostar forma- 
tion. On the other hand, when the molecular cloud has an angular momentum, 
the first core does not disappear even after the protostar formation, and directly 
evolves into the circumstellar disk with a Keplerian rotation. There are two paths 
for the formation of the circumstellar disk. When the initial cloud has a consid- 
erably small rotational energy, two nested disks appear just after the protostar 
formation. During the early main accretion phase, the inner disk increases its 
size and merges with the outer disk (i.e. first core) to form a single circumstellar 
disk with a Keplerian rotation. On the other hand, when the molecular cloud 
has a rotational energy comparable to observations, a single Keplerian disk that 
corresponds to the first core already exists prior to the protostar formation. In 
such a cloud, the first core density gradually increases, maintaining the Keple- 
rian rotation and forms the protostar inside it. Thus, the protostar is born in the 
Keplerian disk. In other words, a massive disk already exists before the protostar 
formation. In each case, the protostar at its formation is already surrounded by 
a massive circumstellar disk. The circumstellar disk is about 10-100 times more 
massive than the protostar in the main accretion disk. Such disks are favourable 
sites for the formation of binary companions and gas-giant planets. 

Subject headings: accretion, accretion disks: ISM: clouds — stars: formation — 
stars: low-mass, brown dwarfs: planetary systems: protoplanetary disks 
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Introduction 



A star is born with a circumstellar disk in the molecular cloud core. Molecular clouds 



have an angular momentum ( Arquilla fc Goldsmithlll986l : iGoodman et al. lllQQSt ICaselli et al 
2OO2I ) that can form a disk in the star-formation process. Thus, disk formation is a natural 
consequence of the conservation of angular momentum in the collapsing cloud core. Ob- 
servations have also supported the existence of a circumstellar disk around young stellar 
objects. The circumstellar disks have been identifie d by spectral energ y distributions, and 
are resolved in their scattered light (see review of IWatson et al.l 120071 ). Thus, there is no 
room for doubt about the emergence and presence of circumstellar disks in the star-formation 
process. 

The circumstellar disk plays important roles in star and planet formation. The planets 
are formed in the circumstellar (or protoplanetary) disks. Recent theoretical studies state 
that planets are formed according to the gravitational instability scena rio with massive disk s 
and according to the core accretion scenario with less massive disks ( iDurisen et al.l 120071 ). 
Thus, the mass of the circumstellar disk determines the planet-formation mode and prop- 
erties of planets. In addition, the protostellar outflows are driven by the circumstellar disk 
( jPudritz et al.l 120071 ). It is considered that these flows are closely rela ted to the angular 
momentum transfer, and they determine the star-formation efficiency (IMatzner fc McKee 
20001). Therefore, it is necessary to investigate the formation and properties of the circum- 
stellar disk to understand both star and planet formation. Nevertheless, so far, the formation 
process of the circumstellar disk in the collapsing cloud core has rarely been investigated. 

In a classical star-formation scenario, it is envisaged that the circumstellar disk begins 
to form after the protostar formation and increases its size with time. This picture is evoked 
by the idea that a small part of the natal cloud having the lowest angular momentum first 
collapses and forms the protostar, while the remainder of the cloud, which has larger angular 
momentum collapses later, and forms the circumstellar disk. However, in this picture, the 
thermal evolution of collapsing (or infalling) gas is not taken into account. Recent star 
formation scenarios including thermal evolution, have shown t hat the transient object that i s 
called the first core appears prior to the protostar formation (IMasunaga &: Inutsuka II2OOOI ). 
In the collapsing cloud core, the gas collapses isothermally for n ^ 10^°cm~^, while it 
collapses adiabatically for n ^ 10^^ cm~^ and forms the first (adiabatic) core at n 10^^ — 
10^^ cm~^. Although the first core is in a nearly hydrostatic equilibrium state, its central 
density gradually incre ases because the first core increases its mass by the gas accretion 



( iMasunaga et al. Ill998l ). When the central density exceeds n ^ 10^^ cm ^ (or temperature 



T ^ 2000 K), the molecular hydrogen begins to dissociate and the central region rapidly 
collapses again (the so-called second collapse). Finally, the collapsing gas becomes adiabatic 
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again and the protostar forms at n ^ 10 cm (lLarsodll969l : iMasunaga &: Inutsuka Il2000l ). 



Because the second collapse occurs only in a small part of the first core, the first core does 
not disappear just after the protostar formation. In a spherically symmetric calculation, the 
first core f alls onto the central protosta r and disappears several years after the protostar 
formation (IMasunaga fc Inutsuka Il2000l ). However, the rotational effect is not included in 
the spherically symmetric calculation. In reality, because the first core is supported by not 
only the thermal pressure but also the rotation, it does not disappear in such a short duration 
without effective angular momentum transfer (jSaigo fc Tomisakall2006l ). 



Recently, iMachida et al.l ( 120101 ) an d iBatd (I20101) pointed out that the first core is the 
origin of the circumstellar disk (see also lBatdll998l ). indicating that the disk appears before 
the protosta r formation and is rn ore massive than the protostar at the protostar forma- 
tion epoch ( Inutsuka et al. 201ol ). Such massive disks tend to exhibit the gravitational 
instability that contributes to the angular momentum transfer through the nonaxisymetric 
gravitational torque, and binary or planet formation through fragmentation. To investigate 
the formation and evolution of the circumstellar disk in the collapsing cloud core, we need 
a three-dimensional simulation that covers both the circumstellar disk and its natal cloud 
with adequate spatial resolution. However, because such a calculation requires a consid- 
erably higher spatial resolution, only a few studies have investigated the formation of the 
circumstellar disk from the prestellar stage. In many other calculations, the evolution of 
circumstellar disks after the main accretion phase was investigat ed with artificial set tings in 
which the disk properties are artificially modeled (see review of iDurisen et al.l 120071 ). How- 
ever, the evolution of the circumstellar disk even after the main accretion phase should be 
investigated from the prestellar stage to give adequate initial condition. Otherwise, artificial 
initial settings may determine the final fate of the circumstellar disk. 

Walch et al.l J2009al Jbl) and iMachida et al.l (120101 ) calculated the formation of the cir- 
cumstellar disk from the prestellar stage (i.e. from the starless molecular cloud core stage) 
up to the main accretion stage in an unmagnetized cloud. However, they did not r esolve 
the protostar itself to perform a long-term calculation of the circumstellar disk. iBatd (Il998[ ) 
investigated the protostar formation with a sufficient spatial resolution up to the protostar 
formation. The cloud evolution in both phases (before and after the protostar f ormat i on) re - 
solving the protostellar radius in an unmagnetized cloud was only calculated by iBatd (120101 ). 
However, that study did not explore the evolution of the circumstellar disk just after the 
protostar formation. 

In summary, so far, the evolution of the circumstellar disk in the proximity of the pro- 
tostar just after the protostar formation has not been investigated with a multidimensional 
simulation resolving the protostar (or protostellar size). If the first core that is formed be- 
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fore the protostar formation falls onto the protostar and disappears for a short duration, the 
circumstellar disk gradually forms a fter the prot ostar formation, as e n visage d in the classical 
picture. Instead, as pointed out by iBatd (I2OIOI ) and iMachida et al.l (120 lOl ). if the first core 
remains longer, it becomes the circumstellar disk that is considerably more massive than the 
protostar. In this study, we calculate the cloud evolution from the prestellar stage up to 
the period after the protostar formation resolving the protostellar radius of < 0.01 AU and 
investigate the formation process of the circumstellar disk in the proximity of the protostar. 
The structure of this paper is as follows. The framework of our models and the numerical 
method are given in §2. The numerical results are presented in §3. We discuss the further 
evolution of the disk in 54, and summarize our results in 55. 



2. Model Settings 



To study the evolution of rotating gas clouds and the formation of circumstellar disks, 
we solve the equations of hydrodynamics including self-gravity: 



0, 



dv 



p{v ■ V)v = -VP 



pV0, 



(1) 

(2) 
(3) 



where p, v, P, and denote the density, velocity, pressure and gravitational poten t ial, re - 
spectively. To mimic the temperature evolution calculated bvlMasunaga fc Inutsuka I (l2000l ). 
we adopt the piece- wise polytropic equation (see, eq [5] of iMachida et al.l 120071 ). In reality, 
the thermal evolution in the collapsing clo ud may have to be investigated with radiation hy- 
drodynamics code. However, iBatd (I2OIOI ) investigated the evolution of the collapsing cloud 
using radiation smoothed particle hydrodynamics code, and commented that the cloud evo- 
lution with radiation hydrodynamics is almost identical to the evolution obtained from the 
calculation with th e barotropic equation of state before the protostar formation (see also, 
Tomida et al.ll2010l ). In this study, we investigate the cloud evolution with the barotropic 
equation of state up to the epoch just after the protostar formation. We think that, in this 
period, the cloud evolution and formation of the circumstellar disk can be safely investi- 
gated by using the barotropic approximation. This treatment reduces CPU time and makes 
it possible to study cloud evolution with sufficiently high spatial resolution. 

As the initial state, we take a spherical cloud with critical Bonnor-Ebert (BE) density 
profile, in which a uniform density is adopted outside the sphere (r > Rc). For the BE 
density profile, we adopt the central density of ric = 10^ cm~^ and isothermal temperature 
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of T = 10 K. For these parameters, the critical BE radius is -Rc = 4.8 x 10^ AU. To pro- 
mote the contraction, we increase the density by a factor of /=1.68, where / is the density 
enhancement factor that represents the stabihty of the initial cloud. With / = 1.68, the 
initial cloud has (negative) gravitational energy twice that of thermal energy. The mass 
within r < i?c is M = 0.8 Mq. Initially, the cloud rotates rigidly with an angular velocity VLq 
around the 2;-axis. We parameterized the ratio of the rotational to the gravitational energy 
(/3o = Er/ Eg) inside the initial cloud. We calculated four different models with different /So- 
Model names, initial angular velocities VLq, and /3o are summarized in Table 1. 

The time step for the calculation becomes increasingly short as the cloud collapses (or 
the central density increases). When the time step becomes extremely short, we cannot 
calculate the evolution of the cloud and circumstellar disk for long durations. To avoid this, 
we adopt a sink at the centre of the cloud after the protostar formation. To model the 
protostar, we adopt a sink around the center of the computational domain when the number 
density n exceeds > 3 x 10^^ cm~^ at the cloud centre. In the region r < rsink = 2 i?©, the gas 
having a number density of n > 3 x 10^^ cm~'^ is removed from the computational domain an d 



added to the protostar as grav ity in each timestep (for details, see lMachida et al.ll2009l . l2010l ). 



Masunaga fc Inutsuka I (120001 ) showed that the protostar forms at n ~ 10^° — 10^^ cm~^ with 
a radius of ~ 2 Rq in the collapsing cloud. This treatment of the sink makes it possible to 
calculate the evolution of the collapsing cloud and circumstellar disk for a longer duration. 
However, in our calculations, the central density never exceeds > 3 x 10^^ cm~^ except for 
the model with /3o = (i.e., non-rotating model). Thus, we calculated the cloud evolution 
up to the end of the calculation without a sink for the model with /3o 7^ 0, while with a sink 
for model with /3o = 0. In each model, we stop the calculation about 4yr after the protostar 
formation. 



T o calcul ate over a large spatial scale, the nested grid method is adopted (IMachida et al. 



2005al . l2006al ). Each level of a rectangular grid has the same number of cells, 128 x 128 x 16. 
The calculation is first performed with five grid levels {I = 1-5). The box size of the coarsest 
grid / = 1 is chosen to be 2^Rc. Thus, a grid of / = 1 has a box size of 1.5 x 10^ AU. A 
new finer grid is generated before the Jeans condition is violated. The maximum level of 
grids is restricted to /max = 21. The / = 21 grid has a box size of 0.14 AU and a cell width 
of 1.1 X 10-3 AU. With this method, we cover eight orders of magnitude in spatial scale. In 
other words, we can spatially cover both the molecular cloud core and protostar itself. 
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Results 



The disk-formation process depends on the rotation rate of the initial cloud, because the 
cloud rotation makes the circumstellar disk in the star-formation process. In this study, we 
calculated the formation of the protostar and circumstellar disk in the molecular cloud core 
for four different models in which the rotational energy of the initial cloud is parameterized 
in the range of < /3o < 10~^. The observations have shown that the molecular clouds 
have rotational energy in the range of 1 0"'^ ^ f^o ^ 0.07 with an average of f3o = 0.02 
( iGoodman et al. Ill993l : ICaselli et al.ll2002[ ). Recent numerical studies showed that the cloud 
with rapid rotation 0r\ >, 10~^ tends to show frag mentation or binary formation before the 
protostar formation (IMatsumoto fc Hanawall2003l ). However, it is difficult to investigate the 
formation process of the circumstellar disk in a binary system. Thus, in this study, we 
limited the rotational parameter in the range of Po < 10~^ to focus on the circumstellar disk 
around the single protostar. In this section, first, we show the cloud evolution for the model 
with a relatively slow rotation rate (/3o < 10^'^). Then, the evolution of a relatively rapidly 
rotating cloud (/3o = 10^'^) is described. The evolution of the cloud with larger rotational 
energy (/3o > 10^^) is discussed in §4.21 



3.1. Evolution of Cloud with a Relatively Slow Rotation 

Figures [1] and [2] show the density distribution and velocity vectors on the z = (Fig. [1]) 
and y = (Fig. [2]) planes for a model with Pq = 10~^ (model 2) with different spatial scales 
tc = 1.2083 yr after the protostar formation. In these figures. Figures [T]a and [2]a cover the 
entire region of the molecular cloud core with a box size of ~ 8000 AU, while Figures [l]f and 
[2]f cover the protostar and a part of the disk with a box size of ~ 0.3 AU. For this model, 
the protostar forms in the collapsing cloud t = 5.064 x 10^ yr after the calculation begins 
(i.e. after the initial cloud begins to collapse). In this paper, we define as the elapsed 
time after the protostar formation where the protostar formation epoch {tc = 0) is defined 
as the epoch when the central density reaches ric = 10^" cm~^ in the collapsing cloud core. 
In addition, we define t as the elapsed time after the calculation begins (or the initial cloud 
begins to collapse). 

Figures [T]a and [2]a show that a molecular cloud with a size of ~ 8000 AU maintains a 
spherical structure even after the protostar formation, because the gas begins to collapse 
in the small area of the cen tre of the clo ud and the outer envelope maintains its nearly 



initial structure, as shown in iLarsoru (119691 ). As seen in Figures [T]& and [2]^, the gas collapses 



spherically towards the centre of the cloud in a large scale of S> 1 AU, because the rotatio n 



hardly affects the dynamical evolution of the cloud at such a scale (IMachida et al. Il2005al ). 
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In the collapsing cloud, the gas becomes optically thick against the dust cooling and the 
first core surrounded by the s hock appears when the number density exceeds n ^ lO^^cm"^ 



(iMasunaga fc Inutsuka Il2000l ). When the cloud has an angular momentum, the first core 
has a disk-like structure. The disk-like structure surrounded by the shock in Figures [He and 
[2]c corresponds to the first core. Figures [T]6-c? show that the gas falls radially towards the 
first core outside the first core, while the gas rotates rapidly inside the first core. This is 
because, inside the first core, the gas collapses very slowly in the radial direction and the 
rotation dominates the radial motion after the gas becomes adiabatic (i.e. after the first core 
formation). 

Inside the first core, there is another structure with a radius of ~ 0.4 AU surrounded 
by the disk-like shock, as seen in panels d and e of Figures [T] and [H This disk-like object 
is formed inside the first core after the protostar formation. In this paper, we call this 
object (i.e., the disk around the protostar) the inner disk. In addition, the protostar with 
a radius of ~ 0.02 AU exists with a nearly spherical structure inside the inner disk, as seen 
in Figures [T]f and [2]f. Thus, after the protostar formation, three different structures appear 
inside the collapsing molecular cloud core: two nested disks (disk-like first core and inner 
disk) and protostar. 

After the first core formation, further rapid collapse is induced in a small central part of 
the first core ow ing to the dis s ociation of molecule hydrogen when the number density exceeds 



n ^ 10^^ cm ^ (jLarsorull969l : iMasunaga &: Inutsuka Il2000l ). In such a collapsing region, the 



gas becomes adiabatic again after the dissociation of molecule hydrogen is completed and a 
protostar with a shock (or the second adiabatic core) appears in the region of n ;^ 10^° cm^^. 
In this model, then, the inner disk appears around the protostar after the protostar formation, 
because the first core has an angular momentum and the rotation of the infalling gas from the 
first core cannot be neglected as it approaches the protostar. Therefore, the disk structure 
is formed inside the first core because of the rotation. Note that the first core appears in the 
collapsing cloud even without rotation, because the first core is supported by not only the 
rotation but also thermal pressure, while no inner disk appears without rotation. Figures [1] 
and [2] clearly show two nested disks inside the molecular cloud core outside the protostar. 

The nested disks can also be seen in Figure [31 in which the first core and inner disk are 
plotted by the orange and red iso-density surfaces, respectively. This figure clearly shows 
that the thin inner disk is surrounded by the thick torus-like disk (i.e., the first core). In the 
figure, the velocity vectors on the equatorial plane are projected onto the bottom wall surface. 
The direction of the velocity vector is suddenly changed at the surface of the first core; the 
gas falls radially towards the centre of the cloud outside the first core, while it falls slowly 
with rapid rotation inside the first core. In addition, inside the inner disk (red iso-density 
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surface), the azimuthal greatly velocity dominates the radial velocity (i.e. v^f, ^ Vr)- Thus, 
on the equatorial plane, the gas inside the inner disk is rotating rapidly without infalling 
towards the protostar. 

Figure H] shows the time sequence of the first core and inner disk after the protostar 
formation, in which each panel has a same box size of ~ 3 AU. In the figure, the outer shock 
surface corresponding to the first core is plotted by the broken red line, while the inner 
shock surface corresponding to the inner disk is plotted by the broken blue line. The first 
core has an ellipsoidal structure at the protostar formation epoch (Fig. HJa). Then, after the 
inner disk appears around the protostar (or after the protostar formation), the first core 
thins in the region just outside the inner disk. In Figure Hf), we can see the sharp drop 
of density distribution in the vertical direction in the range of |x| < 0.5 AU. The thermal 
energy decreases in the centre of the first core, where the equation of state is changed from 
7 = 1.4 to 7 = 1.1 because of the dissociation of molecular hydrogen. Thus, the gas inside 
the first core (or in the central region of the first core) falls vertically onto the protostar or 
the inner disk, because the vertical direction is mainly supported by the thermal pressure 
gradient force not by the centrifugal force. After the protostar formation, the inner disk 
increases its size, while the first core slightly decreases its size, as seen in Figures |lf>-c. 
This is because the region inside the first core gradually collapses to accrete onto the inner 
disk. In addition, the fir st core increases its mass by the gas accretion and shrinks its size 



( ISaigo &: Tomisakall2006l ). Finally, two shocks (or two disks) composed of the shrinking first 
core and expanding inner disk merge to form a single disk (or single shock) = 3.4 yr after 
the protostar formation, as shown in Figure After the merger, the single disk increases 
its size with time. 

To confirm this merger, the radial distribution of the density (a), radial velocity (b), 
azimuthal and Kepler velocities (c) and ratio of the radial velocity to the azimuthal velocity 
(d) are plotted in Figure O in which each quantity is azimuthally averaged. At the protostar 
formation epoch [tc = Oyr; black line), the density profile shows only a single shock at 
r ~ 1.7AU (Fig. |5]a). At the shock surface that corresponds to the surface of the first core, 
the radial velocity suddenly increases to approach zero (Fig. [5]&). Inside the first core (r < 
1.7 AU), the radial velocity gradually decreases again in the range of 0.02 AU ^ r ^ 1.7 AU, 
while it increases in the range of r ^ 0.02 AU. Because the protostar {n ^ 10^° cm~^) has a 
size of ~ 0.01 AU, the infall (or negative radial) velocity approaches zero near the protostellar 
surface. The azimuthal velocity suddenly increases at the surface of the first core (Fig. [5]c). 
Figure EH shows that the radial velocity greatly dominates the azimuthal velocity outside 
the first core, while the azimuthal velocity dominates the radial velocity inside the first core. 
This indicates that the first core is rotating rapidly and collapses slowly. 
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The blue and red lines in Figure |5]a show two nested shocks in the range of 0.1 AU <^r <^ 
2 AU: the outer shock corresponding to the first core is located at r ~ 1 — 2 AU, while the 
inner shock corresponding to the inner disk is located at r ~ 0.1— 0.5 AU. The figure indicates 
that the inner shock is outwardly expanding, while the outer shock is inwardly shrinking. In 
Figure [5]6, we can clearly see two shocks for red and blue lines, in which the radial velocity 
is considerably small except for the very proximity of the protostar (r < 0.01 AU) inside the 
inner shock (or the inner disk). Figure [5]c shows that, inside the first core, the azimuthal 
velocity gradually increases with time and has a nearly Keplerian velocity inside the inner 
disk. In addition. Figure [5](i shows that the azimuthal velocity is over 10 times larger than 
the radial velocity inside the first core. Thus, at these epochs, the gas orbits around the 
protostar with a nearly Keplerian velocity inside the first core. 

Two nested structures (i.e. the first core and the inner disk) merge to form a single 
disk tc = 3.396 yr after the protostar formation (green line in Fig. [5]a). The green lines in 
Figure [5]c indicate that the merged disk is supported by rotation and has a size of ~ 1 AU. 
After the merger, the disk increases its mass and size with time. This disk corresponds to the 
circumstellar disk. At the end of the calculation {tc = 3.5 yr), the circumstellar disk has a 
mass of 0.025 Mq with a size of 1.1 AU, while the protostar has a mass of 3.5 x 10~^ Mq with 
a size of 0.01 AU. Thus, the circumstellar disk is approximately seven times more massive 
than the protostar. 



3.2. Evolution of Cloud with Extremely Slow Rotation and No Rotation 

We also calculated the evolution of the circumstellar disk in the molecular cloud with 
a further lowering of the rotational energy of (Sq = 10~^. However, the formation process of 
the circumstellar disk is almost the same as that of the model with /3o = 10"'^. Also, in this 
model, two nested disks appear just after the protostar formation, and finally they merged 
to form a single circumstellar disk. The circumstellar disk has a mass of 0.016 with a 
size of 0.3 AU, while the protostar has a mass of 4.2 x 10"^ Mq with a size of 0.01 AU at the 
end of the calculation (tc = 3.7yr). 

On the other hand, when the initial cloud has no angular momentum, only a single 
core appears outside the protostar. The time sequence for evolution of the collapsing cloud 
after the protostar formation for the model without rotation {Pq = 0, model 4) is shown 
in Figure [61 The azimuthally averaged distributions of the density and velocity for three 
different epochs corresponding to each panel in Figure [6] are plotted in Figure [71 Also, in 
this model, the first core appears in the collapsing cloud before the protostar formation. 
Figure [Hh indicates that the first core exists without disappearing at the protostar formation 
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epoch. In Figure [71 at the protostar formation epoch {tc = 0, black hne), there are two shock 
surfaces: the outer shock corresponding to the first core is located at r ~ lAU, while the 



inner shock corresponding to the protostar is located at r ~ 0.01 AU. iMasunaga &: Inutsuka 



teOOOl ) showed that the first core falls onto the protostar and disappears for a short duration 
after the protostar formation. Figure [6] indicates that, after the protostar formation, the first 
core shrinks with time and finally disappears in tc ~ 1-3 yr. Figure [7] shows that the shock 
surface corresponding to the first core moves inwardly and merges with the protostar. Thus, 
no circumstellar disk appears for the model without rotation. We also confirmed that, after 
the first core disappears, the density distribution in the infalling envelope is proportional to 
p oc r^^'^ in the range of 0.01 AU < r < 10 AU, while it is proportional to p oc in the 
range of r > 10 AU. This distribution resembles the singular isothermal self-similar solution 
flsiu]ll977h . 



3.3. Evolution of Cloud with a Relatively Rapid Rotation 

As shown in §3.H when the initial cloud has a considerably slow rotation (/3o = 10~^), 
two nested disks (the first core and the inner disk) co-exist for a short duration after the 
protostar formation, in which the inner disk is enclosed by the outer first core. On the other 
hand, when the cloud has a rapid rotation, we cannot distinguish the first core and the inner 
disk in the star formation process. In other words, only a single rotating disk (or rotating 
first core) appears before and after the protostar formation. To investigate the formation of 
the circumstellar disk for the cloud with a relatively rapid rotation, the cloud evolution for 
the model with /3o = 10^^ (model 1) is presented in this subsection. 

The density and velocity distributions 16.7 yr (i.e., tc = — 16.7yr) before the protostar 
formation for model 1 (/3o = 10~^) are shown in Figure ID The figure shows that the first core 
has a size of ~ 8 AU and has already become sufficiently thin before the protostar formation. 
The velocity vectors in this figure indicate that the gas is rapidly rotating inside the first 
core. Figure ini shows the density and velocity distribution for model 1 tc = 1.475 yr after 
the protostar formation with various spatial scales. The entire region of the first core is 
covered in Figure [9]a, while the protostar embedded in the first core is seen in Figure [9](i. At 
this epoch (tc = 1.475 yr), the first core and protostar have sizes of ~ 8AU and ~ 0.01 AU, 
respectively. The arrows in Figures EJa-d indicate that the gas is rapidly rotating inside the 
first core outside the protostar. Although the spiral structure appears inside the first core as 
seen in each upper panel of Figures |9]& and c, no clear shock corresponding to the inner disk 
as seen in Figures [1] and |2] appears in this figure. As shown in Figure [9]fi, the gas falls onto 
the protostar rapidly in the vertical direction and very slowly in the horizontal direction. In 
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addition, the disk-like first core sags downwards in the proximity of the protostar, as seen 
in lower panels of Figures [9]c and d. In Figure [TOl the density and velocity distributions 
before (upper panel) and after (lower panel) the protostar formation are plotted in three 
dimensions. The figure shows that the first core has a sufficiently thin disk-like structure 
before the protostar formation. In addition, just after the protostar formation, the protostar 
is already enclosed by the thin disk. In the bottom wall, we can also see the spiral structure 
that is developed by the excess angular momentum of the disk (or the first core). 

In Figure [TT| the azimuthally averaged density (a), radial velocity (6), azimuthal and 
Keplerian velocities (c) and ratio of the radial to the azimuthal velocity (d) are plotted 
against the radius for three different epochs [the epoch before (black) and after (red) the 
protostar formation and just before the protostar formation epoch [blue)]. In Figure [TTh. 
the shock at ~ 8 AU corresponds to the surface of the first core and remains almost in the 
same position by the end of the calculation. Inside the first core (^ 8AU), the gas density 
increases smoothly towards the centre of the cloud before the protostar formation (black and 
blue lines), while it increases with waves after the protostar formation (red line). This is 
because the spiral structure develops inside the first core after the protostar formation, as 
seen in Figure M Figure HVb shows that the radial velocity suddenly increases to approach 
zero at the surface of the first core (~ 8 AU) in each epoch. Just after the first core forma- 
tion (black line), the radial velocity becomes f ^ ~ inside the first core, indicating that the 
gas collapses very slowly towards the centre of the cloud. As the central density increases, 
the second collapse is induced and the gas collapses rapidly again, as seen in the blue line 
of Figure [TTb . After the protostar formation, the radial velocity inside the first core oscil- 
lates around Vr = 0, because the spiral structure appearing after the protostar formation 
contributes to the angular momentum transfer. Thus, inside the disk (or the first core), a 
fraction of gas can fall onto the protostar, while the remaining gas moves outward with an 
excess angular momentum. However, as seen in Figure [TTl tj. because the azimuthal velocity 
greatly dominates the radial velocity inside the first core, the radial motion is not much 
noticeable inside the first core. 

Figure [Tib shows that, in each epoch, the azimuthal velocity traces a near-Keplerian 
velocity inside the first core. Thus, the first core already has the Keplerian rotation from 
its formation. Then, the first core gradually increases its central density maintaining the 
Keplerian rotation, and the protostar appears in the Keplerian rotating disk (or the Keplerian 
rotating first core). Therefore, as seen in the red line of Figure [Tib, the protostar already 
has a Keplerian rotating disk at its formation. In other words, the protostar is born in the 
Keplerian rotating disk that formed long before the protostar formation. As described in 
§3.1[ also in model 2 (/3o = 10""^), the protostar is formed inside the disk-like structure (or 
the disk-like first core). However, the first core does not reach the Keplerian rotation before 
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the protostar formation (see, Fig. [5]c). Thus, the Keplerian rotating disk (i.e. the inner disk) 
appears inside the first core after the protostar formation and merges with the first core 
to form a single Keplerian rotating disk several years after the protostar formation. This 
difference is caused by the initial rotational rate of the molecular cloud core. Because model 
1 has a larger rotational energy than model 2, the first core has a Keplerian rotation from 
its formation and becomes the circumstellar disk directly after the protostar formation. 

At the end of the calculation {t^ = 1.5 yr), for model 1, the circumstellar disk has a 
mass of 0.062 Mq with a size of 8.5 AU, while the protostar has a mass of 2.4 x 10~^ Mq 
with a size of 0.01 AU. Thus, the circumstellar disk is about 25 times more massive than the 
protostar. The larger ratio of the mass of circumstellar disk to that of the protostar for this 
model is also attributed to the larger initial rotation rate. 



Discussion 



4.1. Further Evolution of the Circumstellar Disk 



The formatio n of the circumstel l ar disk in the collapsing cloud w as investigated in pre 



vious studies fe.g.. lYorke et al.lll993l : lBatelll998l : IWalch et al.ll2009al Jbl: iMachida et al.ll2010f ): 



however, the protostar itself and the region in the proximity of the protostar (r ^ 1 AU) were 
not resolved in such studies. Thus, the formation and evolution of the circumstellar disk near 
the protostar just after the protostar formation have not been understood. In other words, 
such studies did not focus on the actual formation process of the circumstellar disk, because 
the circumstellar disk is smoothly connected to the protostar that has a size of ~ 0.01 AU. In 
this study, we investigated the formation of the circumstellar disk having Keplerian rotation 
from the prestellar molecular cloud core stage through the protostar formation resolving 
the protostar itself, and confirmed that the first core that is formed before the pro tostar 
formation is the origin of the circumstellar disk, as expected in iMachida et al.l (120101 ) . This 
indicates that we need to resolve at least the first core to investigate the formation and 
evolution of the circumstellar disk in the collapsing gas cloud. 

When the molecular cloud core has no angular momentum, the gas inside the first core 
gradually accretes onto the prot ostar and the first core disapp ears in several years after the 
protostar formation, as shown in lMasunaga fc Inutsuka I (120001 ). On the other hand, the first 
core does not disappear and evolves into the circumstellar disk after the protostar formation 
when the molecular cloud core has an angular momentum. 

The protostellar mass Mps, the mass of the circumstellar disk Mdisk, and circumstellar 
disk radius rdisk at the end of the calculation for each model are listed in Tabled! As noted in 
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the table, the mass of the circumstellar disk is much more massive than that of the protostar 
by the end of the calculation (or several years after the protostar formation). The first core 
that is the origin of the circumstellar disk has a mass of ~ 0.01 — 0.1 at its formation, 
while the protostar has a mass of ~ 10"'^ Mq at its formation. The mass of each object (the 
first core and protostar) is determined by the Jeans ma ss at its formation. Th e Jeans mass 
continues to decrease as the cloud collapses (see. Fig. 2 of llnutsuka et al. 1120101 ) . and the first 
core forms at earlier epoch than the protostar. Thus, the first core is about 10 — 100 times 
more massive than the protostar. Then, even in the early main accretion phase after the 
protostar formation, the mass of the circumstellar disk that is directly evolved from the first 
core is much more massive than that of the protostar. 

As shown in §3.1^ the first core evolves into the circumstellar disk in the main accretion 
phase even when the molecular cloud core has a considerably smaller angular momentum of 
~ 10~^ — 10~^ whose value is co r aparable or smaller than the lower limit of observations 
/3o ~ 10"'^ (iGoodman et al. I Il993l : ICaselli et al.ll2002l ). This indicates that, in general, a 
massive Keplerian disk rapidly (or suddenly) appears just after (or before) the star formation. 
In the main accretion phase subsequent to the gas-collapsing phase, both the circumstellar 
disk and the protostar increase their mass by gas accretion. In this study, we could not 
investigate further evolution of the circumstellar disk in the main accretion phase, because 
we calculated its formation and evolution with a suffic iently high-spat i al res olution that has 
a considerably short timestep. On the other hand, iMachida et al.l ( 120101 ) calculated the 
formation and evolution of the circumstellar disk with a relatively coarser spatial resolution 
at the expense of the structure in the proximity of the protostar, and showed that the 
circumstellar disk is more massive or comparable to the protostar by the end of the main 
accretion phase. It is expected that such a massive disk tends to show fragmentation and 
subsequent formation of the binary companions or planet-size objects in the circumstellar 
disk. Even when no fragmentation occurs in the main accretion phase, such a massive disk 
is a favourable site for planet formation after the main accretion ph ase, because gas-giant 
planets can be formed by gravitational instability ( iDurisen et al.ll2007l ). However, we require 
a huge amount of CPU time to investigate the further evolution of the disk after the main 
accretion phase from the molecular cloud core stage. 



4.2. Fragmentation in the Early Collapsing Phase 



In this study, to suppress fragmentation and binary formation, we adopted a relatively 
small initia l rotational energy (0o < 10~^) that is small er than the average of the observations 



~ 0.02 ( iGoodman et al. Ill993l ; ICaselli et al.ll2002[ ). Recent numerical simulations have 
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shown that fragmentation or binary formation occurs just after the fi rst core formation when 
the in itial clou d has a relatively larger rot a tional energy of /3n > 1 0~^ (IMatsumoto &: Hanawa 



20031 . see also iBodenheimer et al. I l2000t iGoodwin et al. 1 120071 ) 



Thus, the mass and size 
of the fragments are comparable to those of the first core (r ~ 1 — 10 AU, M ~ 0.1 — 
0.01 Mq), because the first core fragments to form a few clumps. In addition, although 
the angular momentum is redistributed into the orbital and spin angular momenta after 
fragmentation, each fragment still has sufficient angular mo mentum to form the Keplerian- 
rotating disk inside it, as shown in iMachida et al. I ( l2005bl ). Therefore, it is expected that 
the circumstellar disk forms in each fragment passing through the same evolutional path as 
the non-fragmentation model, and the circumstellar disk in each fragment has a similar size 
and mass to the non-fragmentation models. Note that the disk-formation process may differ 
from non-fragmentation models when fragmentation occurs inside the first core (or in the 
higher- density region), because the properties o f the fragment are different from those of the 
first core. However, also note that, iBatd (119981 ) showed that fragmentation rarely occurs in 
such a high-density region (n ^ 10^°cm~^) in the collapsing cloud core. 

We can expect that, in each fragment, the circumstellar disk with Keplerian rotation ap- 
pears after the protostar formation, and such a disk is more massive than (or comparable to) 
the protostar in the early main accretion phase, as shown in ^ However, as the circumstel- 
lar disks grow, the disks in each fragment may interact. Thus, we do not know whether the 
circumstellar disk in each fragment traces the same evolutional path for non-fragmentation 
model in the main accretion and subsequent phases. To investigate the evolution of the cir- 
cumstellar disk in the binary or multiple systems, further long-term calculations that cover 
binary or multiple systems are necessary. 



4.3. Effects of Magnetic Field 



In this study, we adopted unmagnetized cloud cores as the initial state. However, the 
magnetic field plays a significant role in the star- formation process. For example, frag- 



mentation or binary formati on is suppressed by the magnetic braking (IMachida et al.l 12008 



(IMachida et al.l 



2007 



Hennebelle fc Tevssierll2008f). and the protostellar o utfiows are driven by the magnetic effect 



Hennebelle fc Fromangll2008l ). 



In the collapsing cloud core, however, the magneti c field in the high-de nsity gas region 



n ^ 10 cm ) dissipates by the Ohmic dissipation (iNakano et al. 1120021 ). Thus, a very 



weak field exists inside the first core (or the region around the protostar, see IMachida et al. 
20071 ). In this study, we only calculated the evolution of the circumstellar disk for a short 



duration, in which the circumstellar disk has a density of n ;^ 10 cm , as seen in Fig- 
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ures [H [2] and |9l Therefore, in the early main accretion phase on which this study focuses, 
the magnetic field may hardly affect the evolution and formation of the circumstellar disk, 
because the disk formation occurs in the high-density gas region where the magnetic field is 
decoupled from neutral gas. 

The circumst ellar disk incre a ses it s size and mass by the end of the main accretion 
phase as shown in iMachida et al.l ( l2010l ) , in which the evolution of the circumstellar disk in 
the unmagnetized cloud was investigated. However, when the cloud is s trongly magnetized , 
the magnetic field ma y supp ress further evolution of the disk. Recently, iMellon fc Lil (120091 ) 
and iDuffin fc Pudrita (120091 ) proposed that a larger size disk exceeding ~ 10 AU may not be 
formed in a strongly magnetized cloud core, because the angular momentum of the infalling 
gas is effectively transferred by magnetic braking. As the circumstellar disk grows by gas 
accretion, the low-density gas region {n ^ 10^^ cm^'^) appears, where the magnetic field can 
be coupled with neutral gas and the angular momentum is transferred by magnetic braking. 
In addition, in the magnetized collapsing cloud, the circumstellar disk (or first core) can drive 
the protostellar outflow that also outwardly transfers the angular moraentum from the parent 



cloud, and may suppress the larger disk formation (ITomisakal l2002t iMachida et al. 



2004 



Machida et al.l l2006l . 120091 ). However, in the later gas-accretion phase, the gas with larger 
specific angular momentum falls onto the circumstellar disk, because the outer envelope has a 
larger angular momentum and the gas of the outer envelope falls onto the centre of the cloud 
later. Such infalling gas with larger specific angular momentum may be able to contribute 
to the larger size disk formation, over the magnetic braking catastrophe. To investigate the 
actual size of the circumstellar disk, we need to calculate the disk evolution for the further 
long term in the magnetized cloud core. 



5. Summary 

In this study, we calculated the evolution of the collapsing gas cloud from the molecular 
cloud core until the circumstellar disk with Keplerian rotation appears resolving both molec- 
ular cloud core and protostar itself, and found that the circumstellar disk originated from 
the first core and is much more massive than the protostar in the early main accretion phase. 
In summary, the massive Keplerian disk already exists just after the protostar formation. 
This result is different from the classical picture of the circumstellar disk formation, in which 
the circumstellar disk gradually grows its size and mass after the protostar formation. 

The formation of the circumstellar disk from the prestellar core stage is summarized in 
Figure [121 iii which three different epochs for first core formation, protostar formation and 
main accretion phases are schematically described for three different modes of the protostar 
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and circumstellar disk formation (i.e. for the cloud with rapid [a], slow rotation [b] and no 
rotation [c]). The figure shows that the first core appears before the protostar formation. The 
first core shrinks and disappears in several years after the star formation when the molecular 
cloud has no angular momentum (Fig. [T^). On the other hand, when the molecular cloud 
has an angular momentum, the first core has a disk-like structure at its formation and 
remains without disappearing even after the protostar formation. 

When the cloud has a relatively rapid rotation (Fig. [T^). the first core already has a 
Keplerian rotation from its formation. Then, the first core gradually increases its central 
density, and the protostar appears when the number density exceeds n ^ 10^° cm~^. Thus, at 
the protostar-formation epoch, the protostar is already surrounded by the Keplerian rotating 
disk, which is about 100 times more massive than the protostar. In summary, the Keplerian 
disk formation precedes the protostar formation. In other words, the protostar is born inside 
the Keplerian rotating disk. 

When the molecular cloud has a relatively slow rotation rate (Fig. [T2b). the first core is 
partially supported by the centrifugal force, and mainly by thermal-pressure-gradient force. 
Thus, at the first-core-formation epoch, although the first core has an oblate structure, it 
rotates slowly with a sub-Keplerian velocity. In this case, after the protostar formation, two 
nested cores (or two nested disks) appear as seen in Figure HTb . The outer disk corresponds 
to the remnant of the first core, while the inner disk corresponds to the Keplerian disk 
formed in the proximity of the protostar. In the main accretion phase, the gas of the first 
core gradually falls onto the protostar. Because the sizes of the first core (~ 1 AU) and 
the protostar (~ 0.01 AU) are considerably different, the infalling gas increases its rotational 
velocity as it shrinks in the radial direction according to the angular momentum conservation 
law. As a result, the Keplerian disk (i.e. the inner disk) appears around the protostar. Then, 
the inner disk increases its size, while the first core gradually shrinks. Finally, two nested 
disks merge to form a single disk that has a Keplerian velocity and corresponds to the 
circumstellar disk. At the epoch of the merger, the size and mass of the merged disk are 
almost the same as those of first core at its formation. Note that the first core shrinks before 
merger, while it still has a size of ~ lAU (Fig|5]). Thus, also in a cloud with considerably 
small rotational energy, the massive Keplerian disk already exists just after the protostar 
formation. 

In this study, we showed that a massive disk already exists before the protostar for- 
mation, and the protostar is enclosed by the massive Keplerian circumstellar disk from the 
moment of its birth. Although we only calculated the evolution of the circumstellar disk 
several years after the protostar formation, we can expect that the circumstellar disk in- 
creases its size by gas accretion, as in the classical picture of circumstellar disk formation. In 
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addition, our result indicates that the massive disk (> 0.1 — 0.01 Mq) with a size of > 1 AU 
already exists in the very early phase of the star formation (or before the star formation). 
We expect that such a disk can be detected using future instruments such as ALMA. 

We have greatly benefited from discussions with K. Omukai. Numerical computations 
were carried out on NEC SX-9 at Center for Computational Astrophysics, CfCA, of National 
Astronomical Observatory of Japan, and NEC SX-8 at the Yukawa Institute Computer 
Facility. This work was supported by the Grants-in-Aid from MEXT (20540238, 21740136). 
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Table 1: Model parameters and calculation results 



Model 


/3o 


Oo [s-i 




Mps [Mq] 


Mdisk [Mq] 


''disk lAU] 


1 


10-3 


6.5 X 10" 


14 


2.4 X 10-3 


0.062 


8.5 


2 


10--* 


2.1 X 10" 


14 


3.5 X 10-3 


0.025 


1.1 


3 


10"^ 


6.5 X 10" 


15 


4.2 X 10-3 


0.016 


0.3 


4 










3.4 X 10-3 
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KtALiI ^E*^ 



Fig. 1. — The density distribution [colour and contours) and velocity vectors (arrows) on 
the equatorial plane at the epoch just after the protostar formation [t = 5.064 x 10"^ yr, 
tc = 1.2083 yr) for model 2 {/3o = lO""^) at various spatial scales. Each object (molecular 
cloud core, first core, circumstellar disk and protostar) is identified by an arrow. 
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Fig. 2. — Same as Fig. [T]but on the y = plane. 



Fig. 3. — First core (orange isodensity surface) and inner disk {red isodensity surface) for 
model 2 (/3o = 10~^) at tc = 1.664 yr after the protostar formation are plotted in three 
dimensions. The density distribution on the x = 0, y = and z = plane is projected 
onto each wall surface. The velocity vectors on the z = plane are also projected onto the 
bottom wall surface. 
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(a) tc=0 [yr; 




13 



(b) tc=0-434 [yrl 



20 
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(c) t=1 .664 [yr] 
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Fig. 4. — Time sequence for model 2 [Pq = 10~^). The density distribution [colour and 
contours) and velocity vectors (arrows) on the y = plane are plotted for various epochs. 
In each panel, the inner shock that corresponds to the circumstellar disk is plotted by a blue 
dotted line, while the outer shock that corresponds to the remnant first core is plotted by a 
red dotted line. 




Fig. 5. — Radial distribution of (a) the density, (b) radial velocity, (c) azimuthal velocity and 
(d) the ratio of radial velocity to azimuthal velocity for model 2 (/3o = 10~^) with different 
epochs. The lower-left inset in panel (a) is a close-up view in the range of 0.02 AU < r < 
0.2 AU. The dotted line in panel (c) is the Kepler velocity at each epoch. 

20 18 16 14 12 



(a) tc=0 [yr] (b) tc^Q.795[yr] (c) tc=1 .297 [yr] ^-iOi-k.^i-] 




<^.5 0:) -O.i U/J rj.b -U.i LVM U.i 

X |AU] X [AU] X [AU] 



Fig. 6. — Time sequence for model 4 (/3o = 0) after the protostar formation. The density 
distribution {colour and contours) and velocity vectors (arrows) on the equatorial plane are 
plotted in each panel. The elapsed time tc after the protostar formation is described in each 
panel. 
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Fig. 7. — Distribution of density (upper panel) and radial velocity (lower panel) for model 
4 (/3o = 0). 
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Fig. 8. — Density distribution {colour and contours) and velocity vectors (arrows) on the 
z = (upper panel) and y = (lower panel) planes for model 1 {(3q = 10~^) before the 
protostar formation. 




Fig. 9. — Density distribution {colour and contours) and velocity vectors (arrows) on z = 
(each upper panel) and y = (each lower panel) plane at tc = 1.474 yr for model 1 (/3o = 
10^'^) with different spatial scales. 
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First core (or Circumstellar disk) 




Fig. 10. — Density distribution {isodensity surface) around the centre of the cloud before 
(upper panel) and after (lower panel) the protostar formation. The density distribution on 
the X = 0, y = and z = plane is projected onto each wall surface. The velocity vectors 
on the z = plane are also projected onto the bottom wall surface. 
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Fig. 11. — Radial distribution of (a) density, (b) radial velocity, (c) azimuthal velocity and 
(d) the ratio of radial velocity to azimuthal velocity for model 1 {(3o = 10~^) with different 
epochs. 
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Fig. 12. — Schematic view of the circumstellar disk formation in the collapsing cloud core, 
showing three different evolutional sequences for the cloud with (a) rapid rotation, (b) slow 
rotation, and (c) no rotation. 



